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Introduction. 


The  discovery  of  baroclinic  instability  by  Charney  (1947) 
has  dominated  research  in  meteorology  for  the  last  thirty 
years.  The  large  scale  atmospheric  flow  has  been  regarded  as 
a  baroclinically  unstable  circum-polar  vortex,  the  evolution 
of  which  is  sought  in  the  presence  of  transient  disturbances. 
While  it  has  been  recognized  that  the  mean  flow  of  the 
atmosphere  is  asymmetric  as  a  result  of  the  zonal  distribution 
of  orography  and  heat  sources ,  most  of  the  studies  assumed  the 
mean  flow  to  be  zonally  homogeneous .  General  circulation 
models  took  into  account  the  asymmetric  nature  of  the  boundary 
terms  but  they  were  not  very  successful  in  improving  our  basic 
understanding  of  the  consequences  arising  from  the  asymmetry. 

Yet  such  general  circulation  models  did  indicate  that  topography, 
for  example ,  profoundly  influences  the  behavior  of  the  atmosphere 
CManabe  and  Terpstra,  1974)  as  already  had  been  known 
observationally  for  many  years  C Petterssen ,  1956). 

There  are  two  notable  dynamical  features  which  are  believed 
to  be  related  to  the  presence  of  topographic  and  thermal 
asymmetries.  They  are  the  phenomenon  of  blocking  and  the 
observed  geographical  distribution  of  cyclone  occurrence. 

Blocking  refers  to  persistent  large  scale  flow  anomolies  which 
tend  to  occur  in  certain  geographical  locations.  Such  a 
phenomenon  is  depicted  in  Figure  1  which  shows  that  the  high 
index  westerly  flow  is  blocked  completely  over  the  Atlantic  in 
favor  of  large  meridional  excursions. 

Although  the  blocking  phenomenon  has  been  known  for  many 
years  (Rex  1950  a,b)  only  recently  it  has  become  the  focal  point 
of  intense  research  activity.  Tung  and  Lindzen  (1979) 
suggested  that  atmospheric  blocking  could  be  explained  in  terms 
of  simple  linear  resonance  of  planetary  scale  waves  with 
respect  to  surface  forcing  such  as  continental  elevation  and 
land-sea  differential  heating.  Egger  (1978)  proposed  that 
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blocking  could  be  the  manifestation  of  barotropic  nonlinear 
interaction  among  forced  and  slowly  moving  free  waves.  Charney 
and  Devore  (1979)  suggested,  using  a  highly  truncated  spectral 
model  that  blocking  could  be  one  possible  quasi-stable  equi¬ 
librium  state  of  the  atmosphere.  Kalnay  and  Merkine  (1981) 
performed  numerical  simulations  with  quasi-geostrophic  barotropic 
flows  in  an  open  channel.  The  results  revealed  that  when  the 
system  was  repeatedly  excited  at  some  upstream  location  by 
localized  disturbances  an  effectively  time  independent  response 
occasionally  emerged.  The  non-linear  interaction  of  this  field 
with  a  pre-existing  steady  asymmetric  flow  that  was  generated, 
for  example,  by  localized  topography  or  by  potential  vorticity 
sources  led  to  new  steady  state  configurations  including  block¬ 
ing.  Motivated  by  this  work  Merkine  (1980,  1981)  studied 
analytically  the  phenomenon  of  zero  group  velocity  resonance  of 
Rossby  waves  and  obtained  results  in  agreement  with  the  work  of 
Kalnay  and  Merkine  (1981).  In  particular,  he  showed  that  the 
‘  resonance  phenomenon  can  generate  intense  currents  which 

possessed  closed  circulations  in  agreement  with  the  split  jet 
.  configuration  of  some  blocking  events. 

The  possible  relation  between  the  blocking  phenomenon  and 
synoptic  scale  transients  is  implied  in  the  barotropic  model  of 
Kalnay  and  Merkine  (1981).  Recently  Illari  and  Marshal  (1983) 
have  also  implicated  the  inhomogeneous  eddy  fluxes  arising  from 
synoptic  scale  transients  in  the  maintenance  of  blocking 
patterns.  If  such  connection  exists  then  it  is  important  to 
understand  what  determines  the  geographical  distribution  of 
cyclone  occurence.  But  there  are  other  and  more  obvious  reasons 
for  it,  namely  the  improvement  of  local  forecasts  in  regions  of 
intense  cyclogenetic  activity  such  as  the  western  Mediterranean 


The  problem  of  understanding  local  instabilities  is  quite 
difficult.  The  localization  of  the  phenomenon  requires  dealing 
with  a  spectrum  .of  waves  simultaneously  and  at  the  same  time 
taking  into  consideration  the  fact  that  the  basic  flow  is  zonally 
varying.  However,  if  we  temporarily  ignoie  the  zonal  variability 
of  the  basic  flow  but  accept  the  fact  that  baroclinicity  is  more 
enhanced  in  certain  regions  then  the  concept  of  absolute  insta¬ 
bility  (Merkine  1977,  Merkine  and  Shafranek  1980),  which  is 
discussed  in  more  detail  later,  may  be  applicable.  Indeed, 
Pierrehumbert  (1983)  considered,  in  the  context  of  the  two  layer 
model,  the  simplified  problem  of  the  evolution  of  baroclinic 
instability  with  infinite  meridional  scale  in  a  zonally  varying 
basic  flow  and  showed  that  the  eigenvalue  of  the  local  instability 
is  determined  by  the  constraint  of  absolute  instability  evaluated 
at  the  location  of  maximum  baroclinicity.  The  stability  of  the 
zonally  non-hcmogeneous  atmospheric  fields  of  January  and  July  of 
the  Northern  Hemisphere  was  investigated  numerically  by 
Fredriksen  (1983)  and  considerable  agreement  was  found  between 
the  observed  geographical  distribution  of  the  synoptic  scale  eddy 
heat  flux  and  the  distribution  obtained  from  his  numerical  study. 
However,  basic  understanding  of  the  underlying  physics  is  still 
lacking.  The  various  theoretical  aspects  related  to 
Fredriksen's  work  are  discussed  most  thoroughly  by  Pierrehumbert 
(1983) . 

The  fundamental  difficulty  associated  with  stability 
analysis  of  zonally  non  homogeneous  basic  states  is  the  non¬ 
separability  of  the  zonal  coordinate  in  the  equations  governing 
the  linear  evolution  of  the  perturbation.  However,  the  same 

difficulty  arises  at  an  earlier  stage,  namely  in  the  determination 
of  the  basic  state.  Consequently  only  a  few  basic  state  analyti¬ 
cal  solutions  of  the  non-linear  equations  can  be  found  in  the 
literature.  When  the  forcing  and  the  velocity  field  at  infinity 
are  assumed  to  be  independent  of  the  meridional  direction  the  non¬ 
linear  terms  of  the  quasi-geostrophic  equations  vanish  identically 


and  finite  amplitude  solutions  which  are  independent  of  the 
meridional  direction  can  be  found  by  solving  linear  equations. 
This  property  was  utilized  by  Merkine  (1982)  who  studied  the 
stability  of  such  fields.  The  results  were  novel.  They  showed 
the  existence  of  zero  group  velocity  instability  whereby  an 
unstable  non-propagating  wave  packet  developed.  The  unstable 
wavepacket  pivots  about  long  zero  group-velocity  waves  which 
extract  energy  from  the  non-homogeneous  basic  state  but  are 
unable  to  radiate  it  to  infinity.  The  mechanism  is  similar  to 
resonant  triad  interactions  with  one  forced  component.  The 
instability  is  stronger  and  zonally  more  localized  when  the 
zonal  non-homogeneity  is  stronger  hence  demonstrating  in  a 
unique  way  the  subtle  influence  exerted  by  zonal  non-homogeneous 
effects  on  the  dynamics  of  quasi-geostrophic  systems. 

Deriving  analytically  finite  amplitude  steady  state 
solutions  which  vary  zonally  and  meridionally  is  very  difficult. 
Nevertheless,  when  the  conditions  for  zero  group  velocity 
resonance  are  satisfied  simple  analytical  solutions  can  be  found. 
The  barotropic  rectified  currents  of  Merkine  (1980,  1981)  is  one 
such  example.  The  currents  resemble  typical  atmospheric  flow 
patterns  with  strong  meridional  variability  and  weak  zonal 
dependence.  These  currents  proved  very  useful  for  stability 
studies.  The  investigations  of  Merkine  and  Balgovind  (1983)  and 
Merkine  (1983)  reveal  the  existence  of  unstable  localized  baro¬ 
tropic  wavepackets  whose  spatial  structure  and  eigenfrequencies 
depend  on  two  parameters  which  measure  the  degree  of  supercriti¬ 
cality  (the  strength  of  the  horizontal  shear)  and  the  zonal 
length-scale  of  the  shear  region.  The  results  indicate  that  the 
structure  of  the  instability  is  determined  by  conditions  that 
ensure  the  decay  of  the  wavepacket  at  infinity  and  the  transition 
from  long  to  short  waves  across  a  turning  point  (critical  layer) 
region  which  is  controlled  by  non-parallel  effects.  The  control¬ 
ling  influence  exerted  by  the  weak  non-parallel  effects  on  the 
evolution  of  the  instability  demonstrated  the  weakness  of  the 
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basic  state  zonal  homogeneity  assumption  which  can  be  used,  away 
from  critical  layers,  as  a  diagnostic  tool  only. 

The  stability  study  of  Merkine  and  Balgovind  (1983)  dealt 
with  barotropic  fields.  This  is  a  simplifying  assumption.  The 
atmosphere  and  oceans  are  highly  baroclinic  as  a  result  of  the 
differential  solar  heating.  The  consequence  is  two-fold:  the 
basic  state  is  three-dimensional  and  the  energetics  of  the 
instabilities  is  controlled  by  baroclinic  and  barotropic 
processes  which  are  of  comparable  importance.  The  determination 
of  three-dimensional  basic  states  is  naturally  more  complicated  than 
in  the  corresponding  barotropic  problem  and  one  approach  is  to 
study  the  observed  mean  fields  of  the  atmosphere.  This  is  the 
approach  taken  by  Fredriksen  (1983)  who  was  able  to  obtain 
realistic  results.  The  difficulty  with  such  an  approach  is  that 
the  basic  state  is  rather  complicated  and  there  is  no  way  to 
determine  the  influence  of  various  dynamical  factors  on  the 
evolution  of  the  instability. 

The  approach  described  in  this  report  is  different.  We 
extend  Merkine' s  (1980,  1981)  resonance  studies  to  two  layer 
baroclinic  systems  and  derive  analytically  finite  amplitude 
three-dimensional  fields.  The  horizontal  shear,  vertical  shear 
and  zonal  variability  of  these  solutions  are  controlled 
independently.  This  gives  us  much  flexibility  in  assessing 
separately  the  influence  of  barotropic  and  baroclinic  processes 
on  the  evolution  of  trapped  instabilities  in  zonally  varying 
mean  flows.  The  two-layer  model  and  some  of  its  properties  are 
described  in  Section  2.  The  three-dimensional  steady  state 
solutions  are  described  in  Section  3.  The  general  stability 
problem  is  formulated  in  Section  4,  while  Section  5  formulates 
the  linear  stability  problem.  The  necessary  conditions  for 
instability  based  on  the  parallel  flow  assumption  are  stated  in 
Section  6  and  applied  to  our  particular  basic  state.  The  numerical 
model  is  described  in  Section  7  and  some  preliminary  results  are 
presented.  Section  8  outlines  the  anticipated  future  progress  of 
the  research. 
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2. 


The  model. 


We  consider  the  two-layer  model  for  a  slightly  viscous 
quasi- gees trophic  flow  on  a  beta  plane  as  described  by  Pedlosky 
(1970).  Our  particular  system  consists  of  two  layers  of  homo¬ 
geneous,  immiscible  fluid,  with  equal  undisturbed  depths  confined 

vertically  (in  the  z-direction) and  meridionally  (in  the  y- 
direction)  by  rigid  horizontal  boundaries.  The  system  is 
unrestricted  zonally  (in  the  x-direction) .  The  fluid  density 
of  the  upper  layer  is  slightly  less  than  that  of  the  lower  layer 
so  that  the  Boussinesq  approximation  can  be  invoked.  The  system 
rotates  about  the  vertical  axis  with  an  angular  velocity  which 
is  a  linear  function  of  the  meridional  direction  as  implied  by 
the  beta-plane  approximation.  Centrifugal  effects  are  assumed 
negligible  so  that  in  the  absence  of  relative  motion  the  fluid 
interface  is  approximately  level.  Consistent  with  the  quasi- 
geostrophic  formalism  viscous  effects  are  confined  to  thin 
boundary  layers  adjacent  to  the  rigid  boundaries.  The  interface 
is  assumed  inviscid  and  surface  tension  effects  are  ignored. 

We  denote  upper  layer  fields  by  the  index  1  and  lower  layer 
fields  by  the  index  2  and  decompose  the  total  streamf unction  ¥ 
in  the  following  way: 

ft  1  *  *  9  ft 

*1  =  -  U  +  J  e)U  y  +  (LVt )ip1 

(2.1) 

*  „  1  *  *  ?  ft 

v2  =  "  (1  “  7  £)u  y  +  <L  /t)^. 

The  asterisk  denotes  dimensional  variables.  The  first  term  on 
the  R.H.S  of  (2.1)  describes  a  vertically  sheared  uniform  flow. 
The  characteristic  velocity  is  U  and  £  is  the  non-dimensional 
shear  parameter.  The  second  term  on  the  R.H.S  of  (2.1) 
describes  the  contribution  to  the  total  streamfunction  arising 
from  the  presence  of  localized  potential  vorticity  sources  Sn 
such  as  heat  sources,  for  example,  whose  characteristic  time 
scale  is  t.  L  is  the  dimensional  distance  between  the  meridional 
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walls.  In  the  absence  of  sources  the  stationary  flow  field  is 
everywhere  westerly  and  0  <  e  <_  2. 

The  non-dimensional  quasi-geostrophic  equations  governing 
are 

+  (1  +  7  e)  73eH?2*i  +  F<*2  “  ^i5)  +  +  Fe)t^ 

+  r72t|>1  =  -  AJ(^1,V2^1  +  F(iJ>2  -  ^1>)  +  Sx  (2.2) 

+  (1'7  e)  (^2  +  F(*i  “^2^  +  (e'Fe)3l!r 

+  rV2ij;2  =  -  AJ(*2,V24>2  +  F^-^))  +  S2 

* 

where  L  and  U  are  the  reference  scales  for  non-dimensionali- 
zation.  J  is  the  Jacobian  of  two  function  and  is  defined  as 
J(f,g)  =  FxKy“Fygx*  In  addition  to  e  (2.2)  contains  the 
following  non-dimensional  parameters. 

0  =  0'  L2/U* 

r  =  (vfQ)5sL/DU* 

F  =  2L2f2/(gD(p2-Pl)/2 p2) 

A  =  L/U*t 

where  fQ  is  the  Coriolis  parameter  and  0l  its  gradient  at  the 
reference  latitude.  D  is  the  depth  of  the  system,  v  is  the 
kinematic  viscosity,  g  is  the  gravitational  acceleration  and 
Pn  is  the  density.  0  and  F  measure  the  importance  of  the  beta 
effect  and  stratification,  respectively,  r  measures  the  spin- 
down  effect  arising  from  the  secondary  circulation  induced  by 
the  top  and  bottom  Ekman  layers.  The  meridional  walls  serve 
to  restrict  the  lateral  scale  of  the  domain  and  it  is  immateri  1 
whether  they  are  considered  as  viscous  or  inviscid  since  they 
do  not  induce  a  secondary  interior  circulation.  (Greenspan, 
1968).  A  is  the  parameter  of  non-linearity.  It  measures  the 
strength  of  the  potential  vorticity  sources.  When  A  0  the 


potential  vorticity  sources  are  very  weak  and  non-linear 
effects  can  be  ignored. 

It  is  convenient  to  separate  barotropic  and  baroclinic 
effects.  Consequently  we  write 


*»  *  *  *2> 


f(*l  -  *2) 


(2.3) 


where  and  if/^  denote  the  barotropic  ;  and  baro¬ 

clinic  mode s, respectively .  By  adding  and  subtracting  the  two 
equations  of  (2.2)  and  using  (2.3)  we  obtain  the  following  two 
equations 


&  ♦  »  *r  ♦  +  I  «•  4 


-  AJ ( ijj  ,V24)„)  -  AJ (\p  ,V2\p  )  +  S 
m  m  Tt  T  m 


+  '  2FM  *  6  TO1  *  rV2*T  *  7  E  A  V\ 


(2.4) 


+  F  e 


in?* 


m 


■AJ(^m,V2i|)T)  -  AJ(<Pt  ,V2^m) 


+  2FAJ(^m,i pr)  +  St 

where 

Sm=7(SltS2>  -  Sx  “  7(S1  -  S2>  (2-5) 

are  the  respective  barotropic  and  baroclinic  components  of  the 
potential  vorticity  sources.  Equations  (2.4)  were  used  by 
Merkine  (1982)  to  study  the  stability  of  fields  induced  by 
localized  potential  vorticity  sources  which  vary  in  the  zonal 
direction  only.  Our  goal  is  to  study  the  stability  properties 
of  more  general  fields. 

Inspection  of  (2.4)  reveals  that  the  vertical  shear  of  the 
uniform  current  which  is  measured  by  e  and  the  non-linear  inter 
action  terms  which  are  measured  by  A  are  the  two  physical 
entities  which  couple  the  dynamics  of  the  barotropic  and  baro¬ 
clinic  modes  of  motion.  When  non-linear  effects  can  be  ignored 
(2.4)  reduces  to 


&  *  &7\  *  6  *  r7\  *  7  E  sA  *  Sm  . 


1  _  3  p>2. 


&  *  &  (V  V2P"t)  ‘83r,rt\4£i  7\  *  Fe  £  «v  =  sr 


1  .  3  n2. 


(2.6) 


and  when  the  basic  uniform  current  possesses  no  vertical  shear 
(2.6)  reduces  to 

(  3  ^  3  ■Vn2.  .  „  ^m  .  _ o  2 »  _  „ 


(JL  +  -4-)v  ib  +  0  tt-2.  +  rV^to  =  S 
'■at  3xJ  ym  p  Tx  m 

a  a  o  3t()  0 

(Tt  +  ^  '  2F*t)  +  B  +  rV  =  ST. 

The  reason  for  the  nomenclature  to  and  to  is  ^lear  now. 

m  t 

free  modes  of  the  first  equation  of  (2.7)  are  barotropic 
waves  while  those  of  the  second  equation  of  (2.7)  are  ba 
clinic  Rossby  waves.  Since 


(2.7) 


,ssby 


(2.8) 


*1  *  *m  *  *t  •  *2  *  ‘  *t  '  (2-8> 

(2.7)  implies  that  ij^  =  ip2  for  =  0  and  ^  =  -  to2  for  =  0. 
When  small  vertical  shear  is  present  tom  and  toT  are  coupled: 

The  barotropic  mode  is  modified  by  baroclinic  effects  and  the 
baroclinic  mode  is  modified  by  barotropic  effects.  For  baro¬ 
tropic  forcing  S  =  0  and  (2.6)  indicates  that  to  =  0(e)  and 
t  2  T 

that  the  baroclinic  modification  of  to  is  0(e  ).  When  S  =0 

m  a  in 

tom  =  0(e)  and  the  barotropic  modification  of  toT  is  0  (e^).  When 
e  is  no  longer  small  the  coupling  between  tom  and  toT  is  strong 
regardless  of  the  nature  of  the  forcing  and  the  identification 
of  tom  and  toT  with  separate  barotropic  and  baroclinic  dynamical 
responses  is  not  meaningful  anymore.  Nevertheless,  we  still 
find  it  convenient  to  refer  to  tom  and  toT  as  the  barotropic  and 
baroclinic  modes  of  motion  since  they  assume  thoir  pure  form 
when  e  0.  When  c  =  0  and  A  t  0  the  barotropic  and  baroclinic 
modes  interact  non-linearly  (see  (2.4))  but  the  spectrum  can  be 
separated  into  pure  barotropic  and  baroclinic  modes.  We  seek 


We  seek  solutions  to  (2.2)  which  result  from  localized  forcing  and 
initial  conditions.  Consequently  we  assume  that  ij>  decays  to 
zero  at  infinity  implying  that  the  condition  of  vanishing 
meridional  velocity  on  the  meridional  boundaries  can  be  stated 
as 


5  Vi  =  'J'  =  ’J'  =  0 


m 


on 


y  =  0,1. 


(2.9) 


The  dynamical  nature  of  the  response  of  the  system  depends 
crucially  on  the  characteristics  of  the  free  waves.  Hence,  we 
seek  solutions  of  the  form 


_  i 

■  1+Y  l 

i|) 

r  t  ' 

2 

.  1-Y  . 

.  (kx-  tot) 

e  sin  (miry)  +  c.c.;  m  =  1,2,.. 


(2.10) 


where  c.c.  denotes  complex  conjugation  and  substitute  (2.10) 
into  (2.6)  with  =  ST  =  0  to  obtain  the  following  dispersion 
relation  and  vertical  structure  of  the  waves 


c  .  1  .xfe-c-V  a  t  ^(eW/k)2]’*  (2-11) 

K  +2F  K  2K/0<r  +  2F) 


T  *  ^FE|irK?/k  (2.12 

F  (1  +  J£  -  c)F 

where  =  k^+m^ir^.  Equations  (2.11)  -  (2.12)  are  identical 
to  the  corresponding  equations  of  Pedlosky  (1970). 

In  the  atmosphere  and  oceans  the  radius  of  deformation 

It 

=  (gD(p2~P  ^)/(2p  2 ))  v  f  q  is  comparable  to  (U  /(S’;  which  is 
the  typical  length  of  stationary  Rossby  waves  and  consequently 
we  set  F  =  0.  Baroclinic  instability  sets  in  whenever  the 
shear  parameter  exceeds  a  certain  threshold  value.  When  r  =  0 
the  neutral  stability  curve  is  given  by 


c2  *  a&l _ 

K4(»62-K1*) 

and  the  critical  value  of  the  shear  is  f  e [ 
the  neutral  stability  curve  is  given  by 


( >.13) 


1.  When  r  i  0 


2  .  4  [  6  4  K2r2l 

(23-K2)  *-K2(K2+8) 2  k2 


(2.14) 


which  does  not  reduce  to  (2.13)  when  r  -*■  0 .  This  singular 
behavior  was  observed  first  by  Holopainen  (1961)  who  also 
noticed  that  small  r  destabilizes  the  system.  The  threshold 
value  of  e  derived  from  (2.14)  is  equal  to 

o 

|e|  tat  0.908333  +  0(r  )  which  is  less  than  the  above  critical 
value  for  the  strictly  inviscid  system.  In  both  cases,  however, 
the  spectrum  of  the  unstable  waves  is  restricted  to  0  <  K  <26. 
The  waves  to  be  destabilized  first  when  r  =  0  and  r  0  are 
given  by  K2  =  2^6  and  K2  =  (l+3ls)6/2  +  0(r2),  respectively. 

When  the  shear  is  subcritical  and  the  frictional  effects 
secondary  in  importance  each  branch  of  the  dispersion  relation 
resembles  a  cubic.  Both  branches  are  antisymmetric  functions 
of  k  about  k  =  0.  The  number  of  real  roots  of  u  =  0  for  each 
branch  depends  on  the  sign  of  dw/dk  (k  =  0).  When  it  is  posi¬ 
tive  we  have  one  real  root  and  when  it  is  negative  we  have 
three  real  roots.  When  it  vanishes  we  have  the  coalescence  at 
the  origin  of  the  three  real  roots .  This  latter  situation  is 
of  particular  importance.  When  it  occurs  weak  stationary 
forcing  can  trigger  a  large  stationary  response,  dw/dk  is  the 
group  velocity  of  the  waves  and  when  it  vanishes  at  the  origin 
the  energy  pumped  by  the  sources  into  the  long  waves  does  not 
escape  to  infinity  and  a  strong  response  builds:  up.  This  is  the 
phenomenon  of  zero  group  velocity  resonance  discussed  by  Merkine 
(1980,  1981)  when  he  considered  the  response  of  a  pure  barotropic 
system. 
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3.  Steady  state  solutions. 


We  return  now  to  (2.4)  and  iaek  steady  state  solutions 
for  time  independent  localized  sources.  The  problem  is  non¬ 
linear  and  our  approach  is  to  assume  that  non-linear  effects 
are  small  and  to  expand  the  field  variables  in  the  parameter 
of  non-linearity  A.  Consequently  we  write 


ll>  =  ij;0  +  Ail)  *  + .  . 
ym  ym  ym 


\|<T  =  +  a^t  + . . .  , 


(3.1) 


substitute  (3.1)  into  (2.4)  and  equate  to  zero  the'  coeffic¬ 
ients  of  the  various  powers  of  A.  At  the  lowest  order  we 
recover  the  linear  problem  ,  namely 


m' 


r  V2ip° 


vv2rx 


+ 1 


±_ 

3x 


V2*“ 


=  s. 


m 


(3.2) 


*  7 *  - s 


(Recall  that  F  =  6) 

We  follow  Merkine  (1980,  1981)  and  assume  that  the  sources 
have  the  simplified  form 

Sm  =  .u<5(x)  sinmrry  ;  y  =  0,  ±  1 

(3.3) 

ST  =  v6(x)  sinmiry  ;  v  =  0,  ±  1 

where  6(x)  is  Dirac's  delta  function.  This  choice  is  not  as 
restrictive  as  it  seems  to  be .  At  zero  group  velocity  resonance 
which  is  the  dynamically  relevant  case,  a  single  meridional 
wavenumber  dominates  the  response  whose  zonal  structure  consists 
of  very  long  waves.  Thus,  had  we  considered  a  general  source 
function  S  =  f  (x,y)  the  end  result  would  have  been  the  same 


but  the  algebra  .necessarily  more  complicated  (Merkine  1980, 
1981).  y  and  v  need  not  be  equal. 

The  meridional  structure  of  the  forcing  suggests  that 
(3.2)  admits  solutions  of  the  form 

*  'Pm<x)sin(mlTy)  ,  =  ipT(x)sin(imry) .  (3.3) 

The  equations  that  govern  4>m(x)  and  (^(x)  are  suitable 
for  the  application  of  the  Fourier  transform  defined  as 


fCk)  .  — i_  f  P(x)e-iJ<xdx 

(It)*  . 

(3.4) 

oo 

P (x )  =  — i-r-  f  #(k)elkxdk 
(2lO*  J 


where  the  caret  denotes  the  transformed  variable.  It  follows 

that  the  solution  for  $  and  $  is 

m  T 


&  GO  = 
m 


— i-r  (yfk(k2+m2Tr2+8)  +  ir(k2+m2Tr2)l-vke(k24m2Tr2)/2j/D 

(2ir)*  l  L  J 


(3.5) 

$T(k)  =  — |ykeJ^8-(k2+m2iT2)/2j-v^k(B-k2-m2TT2)  +  ir(k2+m2Tr2)  jj/D 

where 

D  =  (l-e2/4)k6  +  {2(l-e2/4)m2ir2  +  e28/2)k4 

+  ((l-e2/4)m47r4  -  8(B-e2m2ir2/2))k2  (3.6) 

-  (r2 2irk)(k2+m27r2)2. 


The  condition  for  zero  group  velocity  of  the  long  waves,  namely 
di)/dk  (k  =  0  ;  r  =  0)  =  0  can  be  determined  from  (2.11).  It  can  be  shown, 
however,  that  it  is  equivalent  to  the  condition 

(1  -  i  e2)m4ir4  -  8(0  -  em2Tr2/2)  =  0  (3.7) 
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obtained  by  equating  to  zero  the  coefficient  of  k  in  the 
third  expression  of  (3.6).  Solving  for  6  we  obtain  that 

6  =  [e2  +  (16  +  e4  -  4e2)Js]m27r2/4,  (3.8) 

which  is  valid  also  for  super-critical  shear  since  the 

2  2 

long  waves  are  always  stable.  When  e  -*•  0  ,  6  =  m  ir  . 

This  is  the  condition  for  zero  group  velocity  of  the  long 
waves  in  the  pure  barotropic  problem  (Merkine,  1980,  1981). 

It  follows  that  (3.8)  is  the  condition  for  zero  group 
velocity  of  the  barotropic  branch  of  the  dispersion  relation 
(2.11).  Equation  (3.7)  admits  a  second  solution  for  6 
corresponding  to  th6  baroclinic  mode  of  th6  dispersion 
relation.'  In  our  case  it  is  negative  implying  that  the 
long  baroclinic  waves  cannot  possess  zero  group  velocity* 

This  conclusion  is  restricted,  however,  to  the  case  8  =  F 
treated  here.  When  $  i  F  the  long  baroclinic  waves  can  have 
zero  group  velocity  for  sufficiently  large  e  (Merkine,  1982). 

The  implication  of  (3.8)  is  that  when  r  =  0  the  inverse 
Fourier  transform  of  (3.5)  does  not  exist  implying  unbounded 
solutions.  This  follows  from  the  fact  that  the  inviscid 
system  operates  at  resonance.  The  energy  pumped  by  the 
source  into  the  longwaves  cannot  be  radiated  to  infinity  or 
dissipated  by  frictional  processes  and  no  stationary  solu- ’ 
tions  are  possible.  Stationary  solutions  exist  when  friction 
is  included  and  we  assume  that  r  <<  1  which  is  applicable 
for  quasi-geostrophic  atmospheric  and  oceanic  systems. 

The  solution  for  4>m  and  <PT  is  determined  by  applying 
Cauchy’ 8  theorem  to  the  inverse  Fourier  transform  of  (3.5). 
The  six  poles  of  the  integrand  are  the  zeros  of  (3.6)  subject 
to  (3.7).  They  are  given  by 


k  =  i  ir  +  o(r) 


<2m4w4)-i/3expi(jr+  ^L)  1/3 

s  — - - - r—r - 9 -  r  + 

'[(2-eV2)ni‘V+e‘^/2] 


o(r1/3)  ;  n  =  0,1,2 


(3.9) 


k  *  *  it  ( 2  -  e2/2)mV  +  e2'8/2]55/(l-e2/4)Js  +  o(l). 
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It  follows  that  the  solution  for  <o  and  <0  is 

m  t 


y(8  +  m2ir2)  -  vem2Tr2/2 


6r5/3mV4 


R(x) 


ye(8  -  m2Tr2/2)-v(8-m27i2) 


R(x)  = 


1/3 

exp(rx/ k^x)  , 


x  <  0 


-2exp(-r1//3k^x/2)cos(~3,5k1r1//3x  +  |it)  ,  x  ^  0 


(3.10) 


k  =  f  2m V  _  11/3 

1  ^  ( 2-e2/2)m2Tr2+e2B/2-‘ 

3  =  [e2+(i6+e4-4e2)5s]m2rr2/4 . 

For  given  e  and  m  the  resonance  condition  (3.8),  which  is  also 
the  last  expression  of  (3.10)  determines  8  and  the  complete 
structure  of  the  solution  is  known. 

The  dimensionless  form  of  the  total  streamfunction  (2.1) 


*1  =  -  (l+e/2)y  +  Ai/^ 
¥2  =  -  (l-e/2)y  +  Ai|>2 


(3.11) 


and  with  the  aid  of  (2.8),  (3.1),  (3.3)  and  (3.10)  it  can  be 
written  as 

fl+e/21  ,  fu[8(l+s)+m2Tr2(l-c/2)]-vE8Mn2ir2(l-e/2)]' 

x  A  *1 

y-“^'7TT 

r  bm  7T  r)  o  o  o 

l-e/2  y[8(l-e)+mV(l+e/2)]+vE6-mV(l+e/2)] 

(3.12) 


R(x)sirn*y  +  0(r  +  A)  >  . 


The  error  estimate  in  (3.12)  consists  of  two  terms.  The  first 
one  reflects  inaccuracies  due  to  the  asymptotic  estimate  of 
the  poles  of  the  linear  solution,  i.e.,  eq.  (3.9)*.  The  second 
term  estimates  the  relative  error  due  to  nonlinear  effects. 

The  interesting  situation  occurs  when  A  =  0(r  ).  Then,  the 

field  induced  by  the  sources  becomes  comparable  to  the  zonal 
flow  and  the  relative  error  estimate  is  0(r  )  and  hence 

small  in  an  asymptotic  sense.  We  have  reached  an  important 
result.  The  solution  of  the  linear  equation  (3.2)  is  in  fact 
a  very  accurate  finite  amplitude  solution  of  the  full  non-linear 
problem  (2.4).  This  surprising  result  which  rarely  happens  is 
a  consequence  of  the  fact  that  we  operate  at  resonance  condi¬ 
tions.  In  the  presence  of  very  weak  forcing  and  dissipation 
the  resonating  solution  is  almost  an  eigenfucntion  of  the  linear 
equations.  In  other  words,  it  is  a  Rossby  wave  which  is  main¬ 
tained  by  weak  forcing  against  weak  dissipation.  But  it  is  well 
known  that  a  single  finite  amplitude  Rossby  wave  satisfies 
exactly  the  non-linear  equation  and  this  explains  the  small 
relative  error  associated  with  the  contribution  of  the  non¬ 
linear  terms  to  the  solution  (3.12). 

The  solution  given  by  (3.12)  describes  a  current  possessing 
vertical  and  horizontal  shear  which  varies  slowly  in  the  zonal 
direction.  The  length  scale  of  the  zonal  variation  is  0(r  ) 

and  this  non-parallel  effect  is  accompanied  by  a  small  meridional 
velocity  which  is  of  0(A/r2/3)  or  OCr^^when  A  =  (Hr2/3).  The 
varying  part  of  the  solution  decays  to  zero  at  large  distances 
from  the  region  of  the  sources .  The  meridional  structure  of  the 
current  is  very  simple  and  it  does  not  vary  in  the  zonal 

♦ 

Note  that  v,hen  the  sources  are  of  general  type  0(A)  is  the 
proper  estimate  of  all  contributions  which  are  not  circttly 
related  to  resonance. 
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2  2 

direction.  When  e  =  0  ,  0  =  m  ir  ,  the  resonance  is  pure  baro- 

tropic  and  regardless  of  the  value  of  v  (3.12)  reduces  to  the 

solution  derived  by  Merkine  (1980)  for  the  barotropic  case. 

From  (3.12)  it  follows  that  0  is  a  monotone  function  of  e.  It 
.  2  2  2  2  2  2 

xs  equal  to  a  u  ,  1.15  m  it  and  2  a  ir  for  e  =  0,  1  and  2, 

ft 

respectively.  For  typical  mid-latitude  values  of  U  «10  m/sec, 

g 

e  1,  L  5.10  m  and  0*  evaluated  at  45°  the  resonance  condi¬ 
tion  is  satisfied  when  m  =  2.  This  result  is  comforting  since 
the  solution  (3.12)  with  m  =  2  corresponds  to  a  slowly  varying 
jet. 

The  derivation  of  (3.12)  is  based  on  the  properties 
of  the  dispersion  relation  for  long  waves.  Such  waves  are 
always  stable  according  to  the  linear  theory  but  is  the  solution 
(3.12)  meaningful  when  e  is  supercritical  such  that  a  part  of 
the  spectrum  is  baroclinically  unstable?  The  answer  to  this 
question  is  found  in  the  concept  of  absolute  instability  intro¬ 
duced  first  into  geophysical  fluid  dynamics  by  Merkine  (1977). 
The  way  a  medium  responds  to  a  given  forcing  depends  on  whether 
the  linear  instability  is  convective  or  absolute.  When  it  is 
convective  the  baroclinically  unstable  transient  set  up  by  the 
initial  value  problem  leaves  the  domain  of  its  excitation  and 
the  long  time  response  at  a  fixed  spatial  point  attains  the  time 
dependent  characteristics  of  the  forcing.  If  the  forcing  is 
steady  the  response  is  steady.  However,  when  the  system 
supports  absolute  instability  a  part  of  the  unstable  transient 
never  leaves  the  domain  of  its  excitation  and  the  long  time 
response  at  any  given  spatial  location  is  dominated  by  an 
exponentially  growing  solution.  The  growth  rate  is  determined 
by  the  imaginary  value  of  w  corresponding  to  the  condition 
dw/dk  s  0  in  the  complex  k-plane.  The  existence  of  absolute 
instability  in  a  given  situation  depends  on  the  parameters  of 
the  problem.  When  we  apply  the  resonance  parameter  constraint 


to  the  criteria  for  absolute  instability  given  by  Herkine  and 
Shafranek,  (1980)  we  find  that  the  steady  state  solution  is 
physically  realizable  unless  c  «  2. 

The  problem  relating  to  the  stability  properties  of  (3.12) 
is  a  different  issue.  When  certain  conditions  for  instability 
are  satisfied  (3.12)  may  become  unstable  to  perturbations. 

Since  our  steady  state  solution  is  three-dimensional  there  is  a 
rich  interplay  between  baroclinic  and  barotropic  processes  both 
influenced  by  zonal  localization  effects.  The  stability  analy¬ 
sis  of  (3.12)  under  a  variety  of  conditions  is  the  main  thrust 
of  the  study.  We  hope  that  the  results  obtained  will  contribute 
to  the  understanding  of  the  difficult  but  important  problem  of 
the  influence  of  zonal  asymmetric  effects  such  as  topography  and 
heat  sources  on  the  dynamics  of  the  atmosphere  and  oceans. 
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4. 


The  formulation  of  the  stability  problem. 

It  is  convenient  to  start  the  stability  analysis  by 
writing  (2.2)  with  3  =  F  in  the  conservation  form 


s  oii<y .  a  n  r  o  -i 

Tt  +  (1+5e  ■“Ty"^  +  “?3r  3yJ[  V  Ai^1  +  g(AiJ;2-Ai{)1)  +  S(l  +  e)y J 


3A<Ji 


=  -  rV*Ai p  +  AS1 


3  Aif» 


(4.1) 


=  -  rV  A\|>2  +  AS2 


which  states  that  the  rate  of  change  of  the  quasi-geostrophic 

potential  vorticity  of  each  layer  following  the  geostrophic 

motion  of  that  layer  is  determined  by  the  potential  vorticity 

sources  and  the  spin-down  effect  of  the  Ekman  layers. 

We  add  now  a  perturbation  to  the  flow  field  such  that 

Ail;  —Aib  +  6y) 
rn  yn  n 


where  6  measures  the  strength  of  the  perturbation.  When  the 

latter  expression  is  substituted  into  (4.1),  while  taking  into 

2/3 

account  the  fact  that  ip  satisfies  (2.4)  to  0(r  )  (as 

n  1/3 

stated  earlier  the  case  of  interest  is  when  A  =  0 (r  ),' 

we  obtain  the  following  system  of  equations  that  governs  the 

evolution  of  the  perturbation  field  <p 


n 


3tt. 


[£.  U*ie*Ul)£*vi  £  ]q,  *  sji 


3y 


6(. 


34^  3q1  3q 

3x  &y  ”  Ty~  Tx 


-)  ;  0(r 


34^  3tt . 
3y  3x~ 


2/3) 


(4.2) 


■|7+(1.le+U2)|-V2  f- 


1  **2 
Jq2  +  37- 

3q 


*k 


3<p, 


3q2 

w 


3(P2 


W~  7x“‘ 


31T2 


34>2  3tt2 
3y  3y 


2) 


=  0(r2/3). 
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qn,  n  =  1,2  is  the  perturbation's  potential  vorticity  defined 
as 


ql  =  7  ^1  +  &<<p2  "  ^l5 

q2  s  V2<P2  +  ®(<pl  " 


(4.3) 


and  irn,  n  =  1,  2  is  the  potential  vorticity  of  the  basic 
state  (3.12)  defined  as 


3^1  2/3 

ir.  =  0(l+e)y  +  — +  $(\f»9-*,)  +  0(r  '  ; 

1  3y  z  1  1 


(4.4) 


3  2&\p  . 

ir,  =  3(l-e)y  +  +  8(<J>,-i{>9)  +  0(rZ/3). 

1  3y2 


The  geostrophic  velocity  induced  by  the  sources  is  given  by 


U 


n 


3A*n  3A^n 

~W  ’  vn  =  ~nr 


n  =  1,2 


(4.5) 


1/3 

where  Un  =  0(1)  and  Vn  =  0(r  ).  Using  the  properties  of 

the  steady  state  solution  (3.12)  we  can  write  (4.4)  as 


=  8(l+e)y  -  (m27r2  +  8 ) Aij^  +  8A^2  +  0(r2//3) 

7T 2  =  8(l-e)y  -  (m27T2  +  8)A^2  +  PAi/^.  +  0(r223) 
the  gradient  of  which  is  given  by 

Vi^  =  [-  (m27T2  +  S)V1+3V2]i 

+  [6(l+e)  +  (m27r2  +  8)U1-8U2]f  +  0(r2/3) 

7 ir 2  =  [-.  (m2ir2  +  8)V2  +  BV1]$ 

♦  t3(l-e)  +  (m2it2+e)U2-0U13$  +  0(r2/3). 


(4.6) 


(4.7) 
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Equations  C3.2)  and  (3.3)  supplemented  by  the  above 
definitions  of  the  basic  state,  the  boundary  conditions 

<Pn  =  0  on  y  =  0,1  ;  n  =  1 , 2  (4.8) 

and  suitable  initial  conditions  constitute  a  well  posed 

initial  value  problem  which  determines  the  evolution  of  the 

1/  3 

perturbation  correct  to  0(r  )  in  non-parallel  effects. 
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The  linear  stability  problem. 

Our  prime  interest  is  in  the  evolution  of  the  perturbation 
field  under  unstable  conditions.  When  the  amplitude  of  the 
growing  perturbation  is  still  small  compared  to  the  basic 
state,  non-linear- effects  are  of  secondary  importance.  The 
most  logical  and  certainly  the  most  common  starting  point  of 
our  stability  analysis  is  the  linear  stability  problem.  The 
linear  dynamics  is  governed  by  (4.2)  with  6=0.  For  the  sake 
of  completeness  we  summarize  now  the  linear  problem  ignoring 
0(r2,/3)  effects. 


[rt* (i+ie+ui}^+vi  wK 
[A+  (1-5e‘U2)^T+V2  ^r]q2 
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qn  =  V2<pn  +  <-l)n0(«P1-4>2) 


(5.1) 


Vn  ~  0  on  y  =  0 ,1 

t  =  0  «Pn  •  <on(x,y) 

Vtt..  =  [-(m2Tr2  +  B)V1+0V„]i  +  [0(l+e)  +  (m2TT2+B)U.-BUo]$ 

1  12  1  (5.2) 

Vtf2  =  [-(m2ir2+B)V2  +  BV1]i  +  [  BU-eHCmV+B^-Bl^]* 

where  U  and  V  are  derived  in  the  usual  way  from  the 
n  n 

x-dependent  part  of  (3.12). 

Non-parallel  effects  enter  the  linear  stability  problem 

through  the  meridional  velocity  V  and  the  zonal  component  of 

n  1/3 

the  potential  vorticity  gradient  St^/Sx  which  are  of  0(r  ) 

and  through  the  dependence  on  the  zonal  direction  of  the  0(1) 
zonal  velocity. 

In  conventional  stability  studies  the  basic  flow  is 
aeeumed  parallel  and  in  the  zonal  direction.  This  amounts 
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to  ignoring  all  0(r  )  effects  in  (5.1).  The  mathematical 

advantage  of  this  assumption  is  that  the  problem  becomes 
separable  in  x  and  normal  mode  solutions  of  the  form 
f^(y)exp  i(kx-wt)  can  be  sought.  The  faith  in  this  approach 
is  based  on  the  notion  that  non-parallel  effects  can  be 
ignored  and  that  the  physically  realizable  solution  is  the 
one  with  maximum  temporal  growth  rate.  Merkine  and  Balgovind, 
(1983)  who  studied  numerically  the  barotropic  analogue  of 
(5.1),  demonstrated  the  inapplicability  of  the  parallel 
theory.  They  showed  the  strong  controlling  influence  of  the 
weak  non-parallel  effects  on  the  evolution  of  the  instability. 
In  fact  all  the  quantitative  predictions  of  the  parallel 
theory  were  proved  to  be  incorrect.  The  only  virtue  of  the 
parallel  theory  was  in  identifying  the  necessary  conditions 
for  instability.  We  adopt  the  same  approach  here,  and  in  the 
next  section  we  discuss  the  necessary  conditions  for  insta¬ 
bility  derived  from  the  parallel  theory. 

Another  important  point  concerns  critical  layers.  The 
linear  parallel  analysis  leads  to  an  eigenvalue  problem 
determining  the  eigenvalue  and  eigenfunction  corresponding  to 
a  given  wavenumber  k.  In  most  circumstances  the  solutions  of 
the  inviscid  parallel  problem  have  a  singular  behavior  at  the 
critical  layer  where  the  phase  speed  of  the  instability  is 
equal  to  the  basic  zonal  velocity.  The  singularity  is 
unacceptable  physically  and  in  the  past  viscous  or  non-linear 
effects  were  used  to  remove  it.  Merkine  (1983)  demonstrated 
that  the  singularity  can  be  removed  through  the  inclusion  of 
non-parallel  effects.  He  also  showed  that  for  the  range  of 
parameters  typical  of  quasi-geostrophic  flows  the  non-parallel 
critical  layer  is  more  important  than  its  viscous  or  non¬ 
linear  analogues. 


6.  The  necessary  conditions  for  instability  of  the  linear 
parallel  problem. 

When  the  weak  non-parallel  effects  are  ignored  in  (4.1) 
and  normal  mode  solutions  of  the  type  discussed  earlier  are 
assumed  two  necessary  conditions  for  instability  can  be 
derived  (Pedlosky,  1979).  In  our  notation  the  two  conditions 


must  be  somewhere  positive  and  somewhere 
negative  in  the  range  of  y. 

(l  +  i  e(-l)n+^  +  Un)  must  be  somewhere  positive 

in  the  range  of  y. 


(6.1) 


(6.2) 


If  none  of  Sn^/Sy  vanishes  then  they  must  be  of  opposite 
signs.  This  is  the  situation  for  baroclinic  instability 
with  no  horizontal  shear  in  a  two  layer  model.  On  the  other 
hand,  if  =  3ir2/3y  as  in  the  barotropic  problem  then 

both  must  vanish  at  some  point  in  order  to  satisfy  the  first 
condition  for  instability.  In  the  general  case  the  two 
possibilities  can  occur  simultaneously  implying  equal  impor¬ 
tance  of  the  baroclinic  and  barotropic  processes.  This  is 
the  situation  typical  of  our  model.  Condition  (6.2)  has  to 
be  satisfied  at  least  at  one  point  in  one  of  the  layers. 

We  proceed  now  to  apply  conditions  (6.1)  -  (6.2)  to  our 
particular  basic  state.  From  (3.11)  -  (3.12),  (4.5)  and 
(5.2)  it  follows  that 


=  8(l+e)  +  (mV  +  BJUj.  -  BU2 


=  B(1-e)  +  (m2Tf2+B)U2  -  BUX 


(6.5) 


where 


I 


r.  _  A  ^ 

6mV 


u  -  A  ^ 


|li(6(l+e)+in^Tr^(l-e/2)  J-v[8-m^Tr^(l-e/2)]|RCx)cosmiiy 

(6.4) 

|p[ 0(l-e)+m2Tr2(l+e/2)  ]+v[  0-m2Tr2(l+e/2)  ] jR(x)cosnnry  . 


R(x),  k^  and  0  are  defined  in  (3.10).  (6.3)  and  (6.4)  reveal  the 

richness  of  the  stability  problem.  The  potential  vorticity  gradients 
depend  on  e  and  m  which  fix  0  and  and  on  A/r2^3,  y  and  v  in 
addition  to  the  parametric  dependence  on  x.  Our  earlier 
choice  of  m  =  2  does  little  to  simplify  the  analysis.  It 
seems  reasonable  to  classify  the  first  necessary  condition 
for  instability  according  to  the  sources  which  can  be 
barotropic,  baroclinic  or  of  mixed  type. 

6.1  Barotropic  sources,  v  =  0. 

From  (6.3)  and  (6.4)  it  follows,  after  some  tedious 
manipulations,  that 


3irl  1 

— -  =  A,  (l+ie+U1 ) 
3y  1  1  1 


2e(0/m2ir2)2+(0/m2Tr2)  +  l-je 
(l+e)(0/m2TT2)  +  l-iE 


(6.5) 


3Tr2  1 
— -=  A9(l-ie+U0) 
3y  ill 


2  2  -2e(0/m2Tr2)2  +  (0/m2ir2)  +  l+je 


A2  =  m  tt 


(1-e)  (0/m27r2)+l+je 


These  expressions  demonstrate  the  proportionality  of  the 
potential  velocity  gradients  in  each  layer  to  the  corres¬ 
ponding  basic  state's  total  velocity.  The  shear  parameter 
E  is  restricted  to  the  range  0  <  e  <  2  for  which  k^  >  0 . 

A2  is  positive  for  e  <  1  and  negative  for  e  >  1.  It 
follows  that  when  the  sources  are  sufficiently  weak  such 
that  the  total  zonal  velocity  in  each  layer  is  of  one  sign 
the  first  necessary  condition  for  instability  is  satisfied 


25 


*v  .  K 

•r . “ifir  itfni’ii'i Mi .  irmtiiiiMr 


when  e  >  1.  This  is  the  same  condition  required  for 
triggering  inviscid  baroclinic  instability  in  a  zonally 
uniform  flow  (see  discussion  following  (2.13)).  We  expect, 
therefore,  that  in  the  case  of  weak  sources  in  the  sense 
discussed  above,,  the  necessary  energy  for  the  developing 
instability  must  come  from  the  available  potential  energy 
of  the  zonally  uniform  part  of  the  basic  state. 

When  e  <  1,  Aj  is  positive  and  the  first  condition  for 
instability  can  be  satisfied  only  if  the  sources  are 
sufficiently  strong  such  that  the  zonal  velocity  reverses 
its  direction  somewhere  within  the  range  of  y.  When  non- 
parallel  effects  are  considered,  this  flow  reversal  is 
equivalent  to  the  appearance  of  closed  circulatiorB  (Merkine 
and  Belgovind,  1983).  It  can  be  shown  that  closed  circula¬ 
tions  appear  first  in  the  lower  layer.  The  maximum  value 

of  R(x)  in  the  expressions  for  U  is  *3  exp( -tt/6  •  3^)  =*  0 . 64 
.  l  i*  1/  3  n  1 

and  it  is  attained  at  ^3  k^r  x  =  ir/6.  It  follows  that  the 
first  necessary  condition  for  instability  is  satisfied 
whenever 
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which  sets  the  strength  of  the  barotropic  potential  vorticity 
sources  for  given  e  and  m.  The  last  expression  indicates 
that  stronger  shear  requires  weaker  sources  for  satisfying 
the  necessary  condition  for  instability. 

In  the  strong-shear  weak-source  case  discussed  earlier 
the  lower  layer  potential  vorticity  gradient  becomes 
negative  everywhere  when  c  >  1.  In  the  weak-shear  strong- 
source  case  just  discussed  the  potential  vorticity  gradient 
becomes  negative  in  isolated  regions  which  appear  first  in 
the  lower  layer.  This  situation  for  which  the  sign  reversal 
occurs  within  a  given  layer  is  typical  of  two-layer  baro¬ 
tropic  processes.  We  expect  that  in  such  circumstances  the 
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enerqy  source  for  the  instability  will  be  the  available 
kinetic  energy  of  the  basic  state.  When  e  >  1  and  the  basic 
state  possesses  closed  circulations  baroclinic  and  barotropic 
processes  must  be  of  comparable  importance  and  the  developing 
instability  should  be  of  a  mixed  type .  This  is  a  rather  complex 
situation  the  stability  properties  of  which  are  not  very  well 


understood. 


6.2  Baroclinic  sources,  p  =  0  , 

When  the  sources  are  baroclinic  the  potential  vorticity 
gradient  can  be  written  as 


Sir. 


n  2  2 

;  =  m  it 


2(8/m2Tr2)2-3(8/m2ir2)  +  1  -  ie 


- T~ 5 - 1 — 

(8/m  it  )  -  1  +  ~c 


(6.6) 


Sir, 


2_  2 


B(l-e)  -  B2U2  $  B2  =  m  u 


2(8/m2iT2  )  2  -  3(6/m2ir2)  +  1  + 
(8/i/tt2)  -  1  -  |e 


We  observe  that  unlike  the  barctropic  case  the  potential 


vorticity  gradients  are  not  proportional  to  the  basic  zonal 
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velocity.  When  e  =  0,  =  0  regardless  of  the  ratio  A/r 


n 


is  the  resonating  part  of  the 


This  result  is  not  surprising, 
basic  state.  The  zero  group  velocity  resonance  is  associated 
with  the  barotropic  branch  of  the  dispersion  relation  which  in 
the  absence  of  shear  cannot  trigger  baroclinic  response.  The 
response  triggered  by  the  baroclinic  sources  when  e  =  0  is 
0(A).  This  is  a  very  weak  response  and  it  is  not  relevant 
for  our  stability  study  since  it  modifies  the  zonally  uniform 
part  of  the  basic  state  by  an  0(A)  only.  Such  small  contribu¬ 
tions  are  not  included  in  our  description  of  the  basic  state. 
Thus  when  e=G,  3un/3y  =  8  and  the  basic  state  is  stable.  In 
the  absence  of  sources  e  >  1  is  necessary  for  instability. 
This  is  identical  to  the  corresponding  barotropic  case  and 
the  instability  is  purely  baroclinic. 


In  the  case  of  barotropic  sources  the  criterion  e  >  1 
holds  as  long  as  Un  is  sufficiently  weak  such  that  no  closed 
circulation  is  induced  by  the  sources.  The  situation  is 
different  in  the  present  case.  A  weak  U  can  alter  signifi¬ 
cantly  the  stability  criterion  and,  in  fact,  we  expect  it  to 
destabilize  the  flow.  This  is  illustrated  most  easily  for  the 
case  of  e  =  1  for  which  3iT2/3y  =  2m2  ir2^  (B2(  e  =  1)  =  -2m27r2j  . 

We  see  that  no  matter  how  small  Uj  is,  its  meridional  structure 
which  is  proportional  to  cos  (miry)  ensures  the  vanishing  of  the 
lower  layer  potential  vorticity  gradient  at  at  least  one  point. 
In  this  case  the  baroclinic  energy  criterion  is  altered  but  the 
energy  source  of  the  instability  is  likely  to  be  the  available 
kinetic  energy  of  the  basic  state.  However  when  e  >  1  and  Uj 
is  very  small  the  instability  should  be  baroclinic  in  nature. 

In  the  combined  case  of  e  >  1  and  strong  Un  the  instability  is 
expected  to  be  of  mixed  type. 

The  one  to  one  correspondence  between  the  closed  circula¬ 
tions  of  the  basic  state  and  the  closed  contours  of  the 
potential  vorticity  gradients  that  we  observed  in  Section  6  does 
not  exist  in  the  case  of  baroclinic  sources.  However,  in  both 
cases  closed  circulations  appear  first  in  the  lower  layer  where 
the  uniform  part  of  the  zonal  flow  is  weaker.  At  the  end  of 
the  report  we  provide  a  few  figures  which  illustrate  typical 
flow  fields  and  the  associated  potential  vorticity  gradients. 

The  above  discussion  has  demonstrated  the  complexity  of 
the  stability  problem  and  there  is  no  point  for  analyzing  now 
the  stability  criteria  for  the  general  case  when  the  sources 
are  of  a  mixed  type.  Further  advancement  in  our  understanding 
of  the  energetics  of  the  instability  can  only  be  achieved  by 
actually  solving  the  stability  problem. 

In  this  section  and  in  the  previous  one  all  the  information 
was  extracted  by  applying  the  first  necessary  condition  for 
instability  to  our  particular  basic  state.  The  second  necessary 
condition  for  instability  can  always  be  trivially  satisfied  and 
hence  it  is  not  very  informative. 
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The  numerical  model  and  a  few  preliminary  results. 

The  stability  problem  formulated  in  Section  4  must  be 
integrated  numerically.  We  describe  briefly  now  the  numerical 
scheme  and  present  a  few  preliminary  results  pertinent  to  the 
linearized  version  of  the  stability  problem.  The  numerical 
code  was  developed  by  Y.  Bar-Sever  who  is  a  graduate  student 
participating  in  the  research  described  in  this  report.  The 
numerical  method  is  an  extention  to  the  two  layer  model  of  the 
numerical  scheme  described  by  Kalnay  and  Merkine  (1981). 

The  numerical  scheme  developed  for  the  initial  boundary 
value  problem  stated  in  Section  4  conserves  potential  enstrophy. 
A  spatial  staggered  grid  is  used  where  the  prognostic  variable 
qn  (the  potential  vorticity  of  the  perturbation)  is  defined  at 
the  center  of  the  grid  and  tp  (the  streamf unction  of  the  per¬ 
turbation)  is  defined  at  the  corners  of  the  grid.  This  allows  a 
direct  implementation  of  the  meridional  boundary  conditions  at 
y  =  0,1  which  require  the  vanishing  of  the  perturbation  stream- 
function.  The  zonal  extent  of  the  numerical  channel  is 
-20  £  x  20  and  the  sources  are  placed  at  x  =  0.  The  time 
difference  scheme  is  the  N-cycle  (Lorenz,  1971)  with  N  =  4. 

The  potential  vorticity  of  the  perturbation  is  updated  at  the 
numerical  zonal  boundaries  of  the  channel  by  extrapolating  it 
linearly  outward  from  the  interior,  i.e.  in  the  direction  of 
the  group  velocity  of  the  waves  which  are  excited  at  the  central 
portion  of  the  channel.  In  order  to  minimize  spurious  reflexion 
of  waves  from  the  zonal  boundaries  of  the  domain  of  integration 
we  used  the  sponge  layer  method  (Kalnay  and  Merkine,  1981) 
whereby  the  time  scale  for  dissipative  effects  is  decreased 
gradually  toward  the  zonal  boundaries .  The  perturbation  steam- 
function  at  the  zonal  boundaries  is  determined  by  imposing  the 
condition  3  <pn/3  x  =  0  (which  amounts  to  linear  extrapolation) 
in  the  definition  of  the  updated  potential  vorticity  evaluated 
at  the  zonal  boundaries.  This  leads  to  a  pair  of  ordinary 
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differential  equations  which  determine  <pn  at  the  zonal 
boundaries.  The  updated  interior  perturbation  streamfunction 
is  obtained  from  the  updated  interior  potential  vorticity  using 
the  Fishpak  Algorithm  541  for  solving  elliptic  systems  subject 
to  the  boundary  conditions  imposed  on  <Pn  as  stated  above.  This 
algorithm  was  developed  at  NCAR  by  P.  N.  Swarztrauber  (1979). 

The  sponge  layer  device  and  the  procedure  for  evaluating 
qn  and  <Pn  at  the  zonal  boundaries  of  the  domain  of  integration 
are  successful  only  if  the  central  portion  of  the  flow  domain 
is  little  affected  by  end  effects.  The  results  of  Kalnay  and 
Merkine  (1983)  demonstrate  that  this  is  indeed  the  case. 
However,  each  numerical  integration  must  be  scrutinized 
separately. 

The  consistency  of  the  numerical  scheme  was  checked  by 
performing  several  test  runs.  For  example,  numerical  integra¬ 
tions  were  performed  for  the  case  of  super-critical  shear,  no 
sources  and  periodic  zonal  boundary  conditions .  This  is  the 
classical  problem  of  two  layer  baroclinic  instability.  The 
initial  perturbation  was  a  localized  pulse.  Baroclinic  insta¬ 
bility  developed.  The  growth  rate  and  phase  speed  of  the  most 
rapidly  growing  mode  agreed  well  with  the  predictions  of  the 
dispersion  relation  (2.11).  In  other  test  runs  the  results  of 
Merkine  and  Balgovind  (1983)  were  reproduced.  In  all  runs  the 
spatial  grid  size  was  uniform  and  equal  to  0.1  which  provides 
a  reasonable  re^lution  for  the  phenomena  under  investigation. 
This  amounts  to  4000  grid  points  for  each  layer.  For  the  range 
of  parameters  considered  so  far  numerical  stability  was 
ensured  by  choosing  0.05  as  the  time  increment  for  every  4 
cycles  of  the  Lorenz  scheme.  In  many  cases  it  is  necessary  to 
integrate  to  about  20  time  units.  Using  the  Technion's  IBM 
3081D  computer  a  typical  CPU  time  for  such  a  run  is  about 
15  minutes.  It  was  impractical  to  check  convergence  by  decreas¬ 
ing  the  grid  size.  We  rely  on  the  conclusion  of  Merkine  and 
Balgovind  (1983)  which  states  that  only  quantitative  but  no 
qualitative  changes  were  observed  when  a  finer  spatial  grid 
was  used. 
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We  proceed  now  to  describe  some  preliminary  results  per¬ 
tinent  to  the  linear  phase  of  the  evolution  of  the  instability. 
In  particular ,  we  consider  the  influence  of e ,  which  is  the  shear 
parameter  of  the  basic  state  at  infinity,  on  the  evolution  of 
the  instability.  The  results  presented  here  correspond  to  posi¬ 
tive  barotropic  sources  i.e.  p  =  1  and  v  =  0  with  A  =  0.1  and 
r=  0.001.  Figures  2  and  3  depict  the  streamfunctions  of  the 
basic  state  corresponding  to  e  =  0.5  and  e  =  1.5,  respectively. 

The  basic  state  describes  a  blocking  configuration.  The  sources 
induce  a  closed  circulation  at  the  center  of  the  channel.  As  a 
consequence,, the  on-coming  westerly  flow  is  split  into  two 
intense  jets  flowing  along  the  meridional  boundaries  of  the 
domain.  For  the  parameters  chosen  the  closed  circulation  appears 
in  both  layers  but  it  is  larger  in  the  lower  layer  which  is 
characterized  by  a  weaker  flow  at  infinity.  The  two  basic  states 
depicted  in  Figures  2  and  3  are  broadly  similar  but  with  one 
difference.  In  the  absence  of  sources  the  basic  state  is  baro- 
clinically  stable  for  e  =  0.5  and  unstable  for  e  =  1.5.  Based  on 
the  discussion  of  Section  6.1  closed  circulation  is  necessary 
for  instability  when  e  =  0.5  but  not  when  e  =  1.5.  Thus  the 
instability  is  expected  to  be  basically  barotropic  when  e=0.5 
and  of  a  mixed  type  where  e  =  1.5. 

The  results  of  the  numerical  integration  reveal  the  exis¬ 
tence  of  an  unstable  localized  wave  packet  characterized  by  a 
single  complex  eigenfrequency ,  u.  For  e  =  0.5  w =  0 . 84+1 . 88i 
and  for  e  =  1.5  u  =  1 . 00  +  2 . 54i  .  It  follows  that  the  increase 
in  the  baroclinicity  of  the  basic  state  enhances  the  instability. 
However,  without  calculating  the  energy  transfer  mechanisms  we 
cannot  be  certain  whether  the  enhanced  growth  rate  is  due  to  a 
release  of  the  available  potential  energy  of  the  basic  state  or 
that  the  increase  in  e  affects  a  greater  release  of  the  avail¬ 
able  kinetic  energy  of  the  basic  state.  The  energy  calculations 
will  be  reported  elsewhere  as  a  part  of  a  detailed  examination 
of  the  stability  problem.  Figures  4-7  suggest,  however,  that 
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the  instability  is  essentially  barotropic.  In  Figures  4  and 
5  we  see  the  streamf unction  of  the  instability  for  e  =  0.5  and 
it  is  very  similar  to  the  barotropic  instability  of  Merkine 
and  Balgovind  (1983).  The  cells  are  bowed  to  the  east  corres¬ 
ponding  to  a  Reynolds  stress  distribution  which  is  positively 
correlated  with  the  horizontal  shear  of  the  basic  state.  The 
absence  of  phase  changes  in  the  vertical  (Figure  5)  indicates 
that  the  instability  is  barotropic.  The  effect  of  the  vertical 
shear  is  more  noticeable  in  Figures  6  and  7  which  correspond  to 
e?1.5.  The  upper  layer  streamfunction  depicted  in  Figure  6a 
is  similar  to  Figure  4a.  It  indicates  release  of  the  upper  layer 
available  kinetic  energy  of  the  basic  state.  The  situation  is 
different  in  the  lower  layer.  The  cells  shown  in  Figure  6b  are 
bowed  in  a  way  favorable  for  transfer  of  kinetic  energy  from 
the  lower  layer  instability  to  the  lower  layer  basic  state.  From 
Figure  7  it  follows  that  the  lower  layer  field  lags  slightly 
behind  the  upper  layer  field  and  this  suggests  that  potential 
energy  is  transferred  from  the  instability  to  the  basic  state. 

It  must  be  emphasized  that  the  above  considerations  are  based 
on  the  parallel  flow  assumption  and  consequently  must  be 
supported  by  a  more  detailed  study.  However,  if  these  conclu¬ 
sions  hold  then  the  increase  of  baroclinicity  can  occasionally 
enhance  instability  through  barotropic  processes. 

One  other  important  feature  revealed  by  Figures  4-7  is  the 
crisis  region  upstream  of  the  maximum.  This  region  is  charac¬ 
terized  by  a  rapid  transition  from  long  to  short  waves.  As 
stated  by  Merkine  and  Balgovind  (1983)  the  dynamics  of  this 
crisis  region  must  necessarily  be  controlled  by  non-parallel 
effects.  This  concludes  the  brief  discussion  of  some  of  our 
preliminary  results. 


8.  Summary  and  prognosis 

In  the  previous  sections  we  have  reported  on  our  research 
efforts  during  the  first  year  of  support  by  the  AFOSR.  Briefly 
stated,  we  have  developed  an  analytical  three  dimensional  basic 
state  whose  zonal  extent  and  baroclinic  and  barotropic  charac¬ 
teristics  can  be  varied  at  will.  This  basic  state  can  be  one 
prototype  model  for  studying  systematically  the  evolution  of 
local  quasi-geostrophic  instabilities  of  a  general  type. 

Although  the  basic  state  is  rather  simple  the  structure  is  such 
that  the  evolution  of  the  instability  must  be  studied  numeric¬ 
ally.  Consequently,  we  have  developed  a  numerical  model  of 
sufficient  flexibility  for  studying  all  the  aspects  of  the 
instability.  We  have  made  a  few  preliminary  runs  and  the  results 
revealed  the  existence  of  a  localized  instability.  During  the 
second  year  of  research  (March  1,  1984  -  February  28,  1985)  we 
shall  carry  out  a  systematic  study  of  all  aspects  of  the  linear 
phase  of  the  instability.  In  particular,  we  shall  examine  the 
way  the  instability  is  affected  by  the  strength  of  the  sources 
and  their  structure  as  well  as  by  the  zonal  length-scale  of  the 
basic  state  and  the  vertical  shear  at  infinity.  The  third  year 
of  research  (March  1,  1985  -  February  28,  1986)  will  be  devoted 
to  the  non-linear  stability  problem.  The  non-linear  study  deals 
with  the  equilibration  of  the  instability  and  the  distortion  of 
the  basic  state.  These  two  related  processes  are  necessarily 
accompanied  by  wave  radiation  from  the  region  of  origin  of  the 
instability. 
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Figure  captions 

Figure  1:  A  blocking  configuration  over  the  Atlantic. 

Figure  2:  The  streamfunction  of  the  basic  state  induced  by  a 
positive  barotropic  source,  i.e.  p=l,  v=0. 

The  shape  of  the  source  is  given  by  6(x)sin2iry . 

In  this  example  A  =  0.1,  r  =  0.001  and  e  =  0.5. 
a)  upper  layer,  b)  lower  layer. 

Figure  3:  The  same  as  in  Figure  2  but  for  e  =  1.5. 

Figure  4:  The  streamfunction  of  the  instability  corresponding 
to  the  basic  state  depicted  in  Figure  2.  a)  upper 
layer,  b)  lower  layer.  The  amplitude  of  the  insta¬ 
bility  is  arbitrary  since  the  problem  is  homogeneous 
The  cells  propagate  eastward  through  a  stationary 
envelope.  The  eigenfrequency  is  w  =  0.84+1.88i. 

The  maximum  of  both  fields  is  rescaled  to  100.  The 
maximum  of  the  lower  layer  field  is  87%  of  the 
maximum  of  the  upper  layer  field. 

Figure  5:  Cross  sections  at  y  =  1/2  of  the  fields  shown  in 

Figure  4.  Full  line  corresponds  to  the  upper  layer 
and  dashed  line  to  the  lower  layer. 

Figure  6:  The  streamfunction  of  the  instability  corresponding 
to  the  basic  state  depicted  in  Figure  3.  a)  upper 
layer,  b)  lower  layer.  The  amplitude  of  the  insta¬ 
bility  is  arbitrary  since  the  problem  is  homogeneous 
The  cells  propagate  eastward  through  a  stationary 
envelope.  The  eigenfrequency  is  u>  =  1.00+2.54i.  The 
maximum  of  both  fields  is  rescaled  to  100.  The 
maximum  of  the  lower  layer  field  is  53%  of  the 
maximum  of  the  upper  layer  field. 

Figure  7:  Cross  sections  at  y  =  1/2  of  the  fields  shown  in 

Figure  6.  Full  line  corresponds  to  the  upper  layer 
and  dashed  line  to  the  lower  layer. 
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